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Abstract 

Highly nonlinear behaviour of a system of discrete sites on a lattice is 
, observed when a specific feedback loop is introduced into models employ- 

^ . ing quantum cellular automata or their real-valued analogues. 

It is shown that the combination of two operations, i.e. i) enhancement 
of a site's value when fulfilling a feedback condition and ii) normalization 
of the system after each time step, produces relatively short-lived spa- 
I tiotemporal patterns whose mean lifetime can be considered as emergent 

. order parameter of the system. This mean lifetime obeys a scaling law 

' involving a control parameter which tunes the "fault tolerance" of the feed- 

back condition. Thus, within appropriate ranges of the systems variables, 
the dynamical properties can be characterized by a "fractal evolution di- 
mension" (as opposed to a "fractal dimension"). 



1 Introduction 
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■ The study of spatiotemporal patterns in dynamical systems has brought some 

general insights into the mechanisms of pattern generation and evolution. Ap- 
plications exist not only for physical or chemical systems exhibiting chaotic 
behaviour |^ but also for models of biological networks |^ . 

Above all, it has turned out that the systems descriptions via discrete maps 
help to avoid complications arising from the analysis of integro-differential equa- 
tions H, p. One important tool for model building is the use of coupled- map 
lattices, either with global couplings as in ||], i.e. 

x(i + l,j) = (l-C)/Wt,j)) + ^5]/(x(t,z)), j = l,...,7V, (1) 

i=l 
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where t denotes a discrete time step, j the index of elements of a one-dimensional 
array and ^ the coupling constant, or with local couplings permitting only 
nearest-neighbour interaction j7|, i.e. 

x{t + l,j) = {l-C)fix{t,j)) + ^{f{xit,j-l)) + f{x{t,j + l))} j^l,...,N. 

(2) 

The mapping function f{x) is usually chosen as nonlinear function, e.g. f{x) = 
1 — ax^ . Thus the systems dynamics is governed by a diffusing part (tuned by 
Q and a nonlinear transformation (tuned by a). 

Another approach consists of the simulation of the systems' evolutions by 
cellular automata where time, space and the values of the state variables 
(sites) are discrete. As opposed to the coupled-map models, the source of com- 
plexity is the large number of degrees of freedom and not a nonlinear map. 

In this paper we choose a different approach to complex pattern formation: 
we first consider the quantum mechanical analogue of Eq. (|2|) , but with a linear 
mapping function. For the sake of probability conservation, a normalization 
procedure after each time step has to be introduced, which on one hand repre- 
sents a source of nonlinearity. On the other hand, it implies a kind of nonlocal 
information spreading all over the array, similar to the global coupling in (|l|). 
This model was introduced as quantum cellular automaton (QCA) some years 
ago and it was also applied as discretized version of the Schrodinger equation 

In its most elementary form the evolution rule reads as 
^it,j) = ^l={^it-l,j)+t6^it-l,j~l)+td*^it-l,j+l)} , ^,SeC (3) 

with S as coupling constant and J\f as normalization factor yielding 

pit,j):^\^it,j)\\ 
Ef=iP(i,j) = l foraUt. 

The normalized values p{t,j) can be interpreted as probability densities for the 
chosen one-dimensional arrays and are used furtheron for presentation of the 
results. This system has been studied thoroughly . It exhibits the emergence 
of exactly predictable coherent large-scale patterns. In fact, the effect of the 
nonlinearity introduced by normalization is not strong enough to prevent an 
exact long-term prediction of its evolution. 

Therefore, to increase the complexity, we provide the system with a memory 
in which the sites' values of the last few hundred time steps are stored. If a 
specific feedback condition applied to each site is fulfilled for any site, the value 
of that particular site is enhanced. Then the array has to be normalized again. 
Thus the necessary conservation of probability is guaranteed by confining all 
possible (continuous) values of the sites within the domain [0, 1]. Consequently, 
the arising patterns built by a group of nonvanishing sites in the spatiotemporal 
plane can be compared with each other, both when \S\ << 1 and \S\ >> 1. 
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In the following, we shall focus on only one specific temporal feedback op- 
eration It will be shown that the complexity of the systems' behaviours 
increases considerably due to this second source of nonlinearity. We think that 
the arising features of self-organization can serve as canditates to understand, at 
least in principle, some aspects of neural activity, particulary in those domains 
where quantum theory would be necesssary to model neural interactions 113] . 



2 Fractal evolution 

For the purpose of illustration, but without loss of generality, we choose the 
parameters of our model as follows. The systems variables for the "undisturbed" 
one-dimensional QCA (cf. Eq. (^) are the size (i.e. the number of sites), 
the initial values of the state variables \I'(0,j), and the mixing parameter 6 = 
\S\e^'^ , (j) & TZ. They are chosen as 

N = 120 

*(0,40) = 0.1, *(0,60) = 0.9, *(0,100) = 0.3 (5) 
S = Sc-il + i) = 0.02(1 -hi). 

The single feedback operation acts on the time coordinate, such that a kind 
of memory is introduced into the system. This operation becomes effective after 
tmem time steps. Accordingly, the memory of the system consists of an array of 
the size t„iem x A^ where the sites' values are stored and cyclically shifted after 
each time step to keep the most recent values of the array. After having reached 
tmem, at cach time step t > tmem it is checked whether the normalized absolute 
value of each site of the array is similar to its value at the time t — tmem ■ The 
similarity check is done by asking whether the present value lies within a specific 
range of the past value. A "relative interval width" e is introduced as crucial 
control parameter of the model: 

< e for all j, and t > tmem (6) 

where e varies between 0.0 (= 0%) and 1.0 (= 100%). 

If the condition (^) is fulfilled, the value of the site under consideration is 
amplified by setting it to Aamp >> 1 

^ {t, j ^ enlianced — ^ainp • ('^) 

Thereafter, the usual normalization procedure is applied. 

Note that this kind of feedback operation simulates a selection process en- 
hancing the values and, consequently, the further evolutions of those sites which 
fulfill the similarity condition. We expect that because of the simplicity of 
the basic model and the additional feedback operation some features of simple 
neuronal systems could be modelled, at least with respect to their principal 
dynamical structures (cf. |T2|). 



^(^,J) 

^{t- tmem,j) 
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Figure 1: Density plot of a QCA with temporal feedback loop accentuating the 
range of the box values from Liow = 10^® to 10^^. Values larger than 10^^ are 
displayed with the maximal box size. The relative interval width e was chosen 
as e = 4.1%. 



The parameters characterizing the feedback operation are chosen as 

tmem — 200 

6 = 0.041 (=4.1%) (8) 

Aamp — 100. 

To have some background against which the arising patterns can be easily 
discerned a threshold value Liow for the probability densities p{t,j) is intro- 
duced. It both makes it possible to represent the patterns (consisting of sites 
with p{t,j) > Liow) in appropriate plots or diagrams and to define a physical 
lower limit for the operational properties of the system which will be discussed 
in detail below. We choose 

Liou, - 10-^ . (9) 

In Fig. (|l]), the evolution of p{t,j) > Liow from time step t — 800 to 
t = 900 with the parameter sets of (||), (||), and (^ is displayed as density 
plot. The range of values of p between Lio^ = 10^^ and 10~^ is indicated 
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Figure 2: Density plot with the same conditions as in Fig. (|l|), but with e = 
25% . The fragmentation of the spatiotemporal patterns clearly has increased. 

with different box sizes, whereas all values larger than 10~^ are displayed with 
the maximal box size. Due to the temporal feedback operation the continuous 
evolution of the QCA as studied in Ref. @ gets disturbed after tmem- The 
arising fragments, operationally defined by a group of spatially and temporally 
connected sites whose values are larger than the chosen limit L/otu, are rapidly 
distributed all over the plane, irrespective of the patterns' locations due to 
the initial configuration. This feature can be understood in the following way: 
if a site fulfills the feedback condition, its value is strongly amplified relative 
to the neighbour's values. Thus it represents a nucleus for further evolution. 
After a certain time, another site will fulfill the time loop condition and will be 
enhanced, whereas the previous evolution will be damped after normalization. 
This behaviour continues until the damping process due to the appearance of 
new nuclei leads to an extinction of the evolving pattern which now remains as 
a fragment in the density plot. 

This behaviour indicates a certain correlation between e as measure for the 
probability of generating new nuclei on one hand, and the lifetimes r of the 
fragments which depend on the frequency of the appearance of new nuclei on 
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Figure 3: Log-log plot of the mean lifetimes f of the fragments versus the relative 
interval width e. Systems parameters are taken from Eqs. (||),(||), and (H). The 
data are fitted with a power-law function. Thus a "fractal evolution dimension" 
characterizing the pattern formation can be derived. 

the other. This correlation can be seen from the comparison of Fig. ([^) with 
Fig. where e has been enlarged up to e = 25% . The fragmentation of the 
evolving patterns has clearly increased, and the mean lifetime of a single pattern 
has decreased. Therefore, a characterization of this new statistical property of 
the system by taking the mean lifetime of the patterns f = f(e) is justified. 
Now the important role of a threshold value for p becomes clear: operationally, 
f, being a kind of "order parameter" of the system, has to be constrained by 
physical limits. 

By varying e between 0.12% and 90%, the plots were analyzed with regard 
to the mean lifetimes f of the patterns. A "pattern recognition" program was 
written to count and analyze each bounded fragment up to i = 5 x 10"* and 
to calculate f (simply as arithmetic mean) and its mean error a. The latter 
includes also a systematical error, mainly due to ambiguities in choosing well 
defined fragments. With the help of parallel visual checks of the plots it was 
estimated as roughly 5% of f. 

The result can be seen in Fig. (^, where the plotted data are shown in a log- 
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log diagram. As a surprising result, the data produce, at least up to e « 35%, a 
linear function, indicating a kind of effective "fractal dimension" of the patterns 
independent of the relative interval width e. The fit was performed with the 
power-law function 

f = a • e'' (10) 

where 

a 2.00 ±0.03, 5== -0.675 ±0.005 (11) 

with a of 14.9 for 52 degrees of freedom. 

We want to emphasize the difference in our use of b as opposed to the usual 
measure of a fractal dimension of an object. Usually, the fractal dimension 
exhibits an "object's" scale invariant features, whereby the process of its deter- 
mination does not change the intrinsic properties of the "object". In our case, 
the scale invariance of the fragmentation of the evolving patterns is generated 
by the repeated application of different values of the relative interval width. We 
call this dynamical fragmentation process a "fractal evolution" of the system. 
The scale invariant property of the emerging order parameter of the system, 
characterized by the mean lifetimes of the patterns, will be denoted as "fractal 
evolution dimension" (De)- In analogy to the usual algorithm |jl^, one gets 
the "fractal evolution dimension" De as 

De = I -b= 1.675 ± 0.005 . (12) 

Note also that the fragmentation ceases (i.e. the lifetime of the pattern at 
the beginning goes to infinity) as e goes to zero (i.e. when the precision of 
the comparison reaches a maximum). This behaviour is plausible, because the 
probability for the feedback condition getting fulfilled is very small for small 
values of e, so that the evolution of the QCA remains nearly undisturbed. 

The deviation of the data points in Fig. (^ from the linear function for 

e ~ 35% can be understood at least qualitatively by looking at the shapes of the 
relative frequencies for the different lifetimes r belonging to a specific e. In Fig. 
(^, they are displayed for three different values of e, i.e for e = 0.9% , 4.1% , and 
25% , respectively. Note that with the default parameter values as introduced 
above no fragments of length t = 1 can occur: assuming the enhancement of 
some site at location j and time t, the value of p{t,j) will be about 1.0 after 
normalization, whereas the other sites will be damped by a factor of about 
TV « ^amp — 10^- Even if, at the next time step, 10 different sites would fulfill 
the feedback condition and thus become enhanced, the normalization factor 
would be of the order A/" ~ 10 x 10**. Therefore, p(t + l,j) would indeed be 
reduced to about 10~^, but it would still remain above the chosen threshold 
of Lio^ = 10~^. Thus, a minimum length of 2 for each arising fragment is 
guaranteed. 

The distributions are all asymmetric, showing poisson-like "tails" propor- 
tional to e""^, but for large e the cut-off at the lifetimes' value of t = 2 becomes 
effective. Thus, no fragments with t < 2 can contribute, and the mean value of 
those lifetimes will be shifted to a higher value. 
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Figure 4: Distribution function of the patterns' lifetimes obtained for tmem = 
200 and for e = 0.9% , e = 4.1%, and e = 25%, respectively. Characteristic 
asymmetric functions are observed. 



We also report briefly on the variation of the systems parameters . It 
turned out that fractal evolution is practically invariant under i) variation of 
the phase of the coupling parameter S (cf. Eq. (||)), ii) the initial values of the 
state variables, and Hi) the choice of other time intervals (than from tmem to 
t = 5 X 10"') for the pattern analysis. 

The variation of the remaining variables in (||) , (H) , and in most cases 
also yields fractal evolution, but with shifted values of a and b in (|l^). Devi- 
ations from fractal evolution are observed in cases of "transient states" of the 
system according to specific values of the parameters. In all cases where such a 
deviation occurs also the shape of the distribution function for f deviates from 
that described above. 

Finally, even the rule (|^) itself can be changed in an appropriate way with- 
out loosing the property of fractal evolution. For example, one can model a 
transition to a "semi-classical" system by restriction to real-valued parameters 
and couplings, i.e. 

^bitj) = —l={^{t - 1, j) + S[^it - 1, J - 1) + ^{t + 1)]} (13) 
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with ^, 5 £ TZ. Also in this case a fractal evolution can be observed. This is 
opposed to the effects of the restriction to real values in the undisturbed case 
[ p^ , where a "real- valued QCA" always leads to trivial final states. 

3 Conclusions 

An implementation of a specific temporal feedback operation in normalized dy- 
namical systems on a lattice with a linear mapping function provides new and 
surprisingly rich evolutionary properties. Specifically, the temporal behaviour 
of QCA evolution with feedback loops selecting and amplifying particular regu- 
larities was investigated. We have found that these regularities are independent 
of the initial point configurations, of the phase of the coupling parameter, and 
of the choice of the time interval whithin which pattern analysis is performed. 

Most importantly, we have found a strict scaling behaviour of the emerging 
spatiotemporal patterns: on a log-log scale, lifetimes f of the patterns' fragments 
are a strictly linear function of the relative interval width e fed back into the 
loop comparing the amplitudes at different evolution times. In other words, 
there exists a "fractal evolution dimension" De oi the patterns whose value is 
independent of the particular choice of e. As e is the relative interval width, De 
is a quantity characterizing all the emerging patterns irrespective of the absolute 
values of the intervals in the temporal feedback loop. Above all, the observed 
features are not restricted to one specific model. It can be concluded that 
discrete, normalized feedback systems with linear mapping functions generally 
exhibit the property of fractal evolution [ p^ . 

Apart from being a "toy model" to study pattern generation in nonlinear sys- 
tems, the observed features of our model are of particular interest for possible 
applications of QCA like modelling processes in the brain as discussed in Ref. 
[ p!2| . As much of the encoding and decoding of neural signalling is a function of 
the "lifetimes" of neural firing, we have seen that our model is of interest because 
the system generating the characteristical "temporal patterns" does not have to 
operate with a "knife's edge" precision. On the contrary, our model gives an 
example of the unimportance of precision in some pattern recognition tasks: It 
just takes any two values for the average lifetimes f at different relative inter- 
val widths e to determine the "fractal evolution dimension" De typical for the 
feedback loop in the QCA, and thus to obtain a measure of f under arbitrarily 
high precision e (that is, within physical limits, of course). Such features may 
be relevant for real-time analog signal processing (as opposed to digital "all- 
or-nothing" processing with correspondingly high precision requirements) which 
allows for high effective noise levels and fault tolerances of neural computations 

[0. 
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